玩转森林图的n种方法,总有一种适合你!

据我个人的不完全调查，可以实现森林图的有：

Meta分析专用软件RevMan

通用统计软件SAS/Sata/R

专用绘图软件GraphPad Prism

大众情人Excel

对于传统的Meta分析里的森林图，在RevMan里填好数据就自动生成了，市面上汗牛充栋的Meta分析培训班也就专门教这个，在此不再赘述。当然，这种传统森林图用Stata的metan命令和R的Meta和Metafor包也能轻松搞定，不太了解的可以戳这篇文章，在此也不再赘述。

Meta 分析的10个问题：从理论概念到编程实践

SAS其实通过GTL也能做出图文混排的森林图，不过比较复杂，我个人不推荐。但你就是好这口的话，我曾修过一个宏，这篇文章里，有过提及，也不再赘述。

多图预警：如何又快又美地用SAS画各种统计图形？

都不再赘述？那今天述啥？今天给大家展示：

如何用大众情人EXCEL搞定期刊风格的亚组分析森林图！

（三）Excel开挂

比如说，顶级医学期刊NEJM的亚组分析森林图的风格如下：

↑NEJM亚组分析森林图的部分截图

分析：难点就在于如何画出HR及其95%CI，并整合到表格中去。

思路：

做好表格，空出Hazard Ratio(95%CI)栏的图的位置。

做好Hazard Ratio(95%CI)的图。

把Hazard Ratio(95%CI)的图贴到表格中去。

是的，这和把大象放进冰箱的思路是一样的。怎么做？转备好截图的疯狂轰炸！

1. 做好表格，留空位置；准备辅助数据，画图。

↑模板示例

辅助数据包括：

标明行号的id列（倒序，以后当做y轴数据用）

效应值及其95%CI(hr, lhr,uhr)

95%CI上下限距离点估计的距离（ler,uer）

无效线列(ref, 通常为0或1，看效应指标是差值还是比值）

2. 插入一个空的散点图。

↑空散点图

3. 右击空图，添加数据，把点估计的效应HR的散点画出来。

注意HR要作为X轴，ID作为y轴。

↑点估计效应的散点图

4. 单击散点，出现加号，仅勾选坐标轴，设置误差线。

注意：13版可能有所差别。

↑设置误差线

5. 用Uer和Ler列分别作为误差线的上限，下限值。

↑继续设置误差线

6. EXEL默认会自动添加固定的纵向误差线，删除他，并对点，线，坐标轴做些设置。

↑删除纵向误差线

↑设置横坐标轴

↑设置效应点，默认圆点，改为正方形

↑初具雏形

7. 添加无效线。右击图后选择数据，Ref列作为X轴，ID列作为Y轴。

↑添加无效线

但是不对啊？这是无效散点啊。

8. 格式化改点为线。

↑格式化无效线

9. 删除Y轴数字标签。调整大小，宽度，移到备用的表格中。

↑八九打扮，九分模样

等等，发现好像哪里不对？是的，要设置透明图片。

10. 设置图片透明。

↑貌似OK

好像还是不对，怎么还有两条竖线线？

11. 设置无边框图片。

（四）皆大欢喜

经过11步大法，成品如下！基本可以出手了。再截图，再导出矢量图片，等等。后续你爱怎么折腾就怎么折腾吧。

哎，等等，怎么又是个11步啊，难道还要过一个双十一啊。哎，……

↑可以出手了

（五）众口难调

不愿意码代码的童鞋，各位临床大大说，「太好了，EXCEL酷毙了，我要投NEJM！」

不愿动手的码农说，「太TMD的麻烦了吧，有没有简单的代码版的啊？」

其实，其实R的forestplot函数也是可以的。当然，其中也是有各种坑，各种小技巧，在此不再赘述。

R版的三行代码效果如下图：

（六）读者福利

最近事儿比较多，写得懒了，前几篇文章放水太多，今天给的这个干货你打几分？打几分？

关注StatsThinking，后台回复「ExlMeta」 我给你推送Excel模板。以后有图，你直接往模板里填数据就可以了。

时间紧了，写得少了，截图太累了，写完这篇我要休息几天。大家且看且珍惜吧！

版权说明：

感谢支持鼓励，形式包括不限于转发，转载，推荐，赞赏。

微信公众平台的 非商业培训机构公众号，无需授权，可直接转载。

网站转载，不得省略作者信息，包括公众号二维码。

内容转载自公众号

StatsThinking

了解更多返回搜狐，查看更多